Ground state of an impurity in a quasi-two-dimensional Fermi gas 
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We study the ground state of a highly imbalanced Fermi gas under quasi-two-dimensional con- 
finement. We find that for typical experimental conditions, the location of the recently predicted 
polaron-molecule transition is shifted to lower values of the vacuum binding energy due to the inter- 
play between transverse confinement and many-body physics. The energy of the attractive polaron 
is calculated in the 2D-3D crossover and displays a series of cusps before converging towards the 3D 
limit. 



Understanding the many-body properties of two- 
dimensional (2D) cold atomic Fermi gases may provide 
insight into puzzling condensed matter systems, such as 
exotic pairing in heavy fermion materials, quasi- 2D or- 
ganic superconductors, and thin film superconductors [1- 
3]. What these systems have in common is that they 
display intriguing behavior when pairing of electrons in 
opposite spin states is frustrated by strong Zeeman mag- 
netic fields, effectively causing an imbalance between 
spin-t and ^ fermions. In this Letter we study the analo- 
gous problem of the highly polarized limit of a quasi-2D 
Fermi gas. 

Recently, several experimental groups have realized 
quasi-2D Fermi gases with ultra-cold atoms [4-10]. In 
these experiments a three-dimensional (3D) Fermi gas is 
confined into a single layer (or a stack of layers) by a tight 
transverse harmonic oscillator potential V{z) = ^mulz^ . 
If the temperature T and Fermi energy ep = h'^kp/2m^ 
where kp is the Fermi wave vector and m is the atomic 
mass, satisfy ksT <C ep <C huz the gas is quantum de- 
generate and collisions can be considered to be quasi-2D 
as transverse degrees of motion are frozen out [11]. 

We consider the highly population imbalanced limit of 
a single spin-^ impurity in a Fermi sea of spin-^ atoms. 
In this limit, it was predicted that the Fermi sea can de- 
stroy pairing (if interactions are sufficiently weak) and 
give rise to a Fermi polaron [12, 13], an impurity dressed 
by excitations of the Fermi sea, in addition to the molec- 
ular state appearing at stronger coupling. From a di- 
mensional analysis of the quasi-2D Fermi gas, the loca- 
tion of the polaron-molecule transition should depend on 
cb/^f and an additional parameter describing the quasi- 
2D confinement (i.e. ep/J^z)- This is to be contrasted 
with the previously studied case of infinitely tight 2D 
trapping [14], whose properties depend solely on the di- 
mensionless interaction parameter cb/^f- As shown in 
Fig. 1, we find that the value of cb/^f where this tran- 
sition occurs strongly depends on the value of ep/fujOz' 
The transition is shifted to lower ^bI^f at higher densi- 
ties, in good agreement with the experimental measure- 
ment of Ref. [10]. The 2D limit considered in Ref. [14] 
is recovered as ep/fujOz 0. The 2D description breaks 
down as the size of molecules (or polarons) approaches 
the harmonic oscillator length of the transverse confine- 
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FIG. 1: (color online). Phase diagram of a highly polarized 
quasi- two dimensional Fermi gas. Below the line, the single 
spin- 4^ impurity exists as a polaron dressed by the Fermi sea of 
spin-t atoms. For large binding energies, the impurity forms 
a molecular bound-state. The location of the transition de- 
pends on the strength of the external confinement: At higher 
density (or weaker confinement) the transition occurs at lower 
tB I^F- The result from the recent experiment by Koschorreck 
et al. [10] is also shown (star). 

ment Iz = ^JhjmuJz' 

We also demonstrate how this quasi- 2D physics is man- 
ifested in the energy of the repulsive polaron of the quasi- 
2D upper branch gas [15, 16]. As we will argue, in the 
repulsive branch, a single effective interaction parameter 
describing quasi-2D collisions provides a sufficient de- 
scription as long as the (non-universal) lifetime of the 
repulsive polaron is large. 

Finally, we present results for the energy of an impu- 
rity in the dimensional crossover regime from quasi-2D 
to 3D, where ^ hhJz- This regime was studied exper- 
imentally in Ref. [6], and also very recently in Ref. [9] 
where features in RF-spectra were interpreted as transi- 
tions between polaronic states. The energy is found to 
exhibit an intriguing behavior, as it has a cusp each time 
the chemical potential of spin-^ particles admits an ex- 
tra harmonic oscillator level. Additionally, the energy is 
found to quickly approach the 3D limit. 
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The low energy scattering of two particles in a two- 
dimensional geometry is described through the 5-wave 
scattering amplitude [17] 
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with q the relative momentum. In the following we work 
in units where = 1. The scattering at low energies is 
controlled by the parameter a2D which has the dimension 
of length. a2D may be extracted in experiments from 
the low energy behavior of the scattering cross section 
o=\!{<l)?l^q. 

Petrov et al. [11] demonstrated how the scattering 
amplitude in a quasi-2D geometry is related to Eq. (1): 
In three dimensions, the properties of an atomic gas in- 
teracting close to a broad Feshbach resonance is charac- 
terized by a single few-body parameter, the scattering 
length ttg. The energy dependence of the quasi-2D scat- 
tering amplitude /oo at energy e = /c^/m < cj^ was found 
to be 
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where incoming and outgoing particles are assumed to be 
in the ground state of the relative motion in the harmonic 
oscillator potential. We use the definition of T [18] 
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As opposed to scattering in 3D, the two-dimensional scat- 
tering always admits a bound state, whose binding en- 
ergy > satisfies Iz/cis = ^{^b/^z)- At low energies, 
|e| <C cj^, the function F takes the form 



^ \n{Bx/n)^^Lx^O{x') (4) 



with B ^ 0.905 [11, 18]. Keeping the first term on the 
r.h.s., the asymptotic low energy expression (1) is recov- 
ered with 



<^2D = /;^\/7r/5exp(- 



(5) 



As demonstrated below, going beyond the leading order 
in the low energy expansion of T{x) is important for 
the understanding of current experiments on Fermi gases 
confined to quasi-2D. 

We now turn the problem of a single | impurity in a 
Fermi sea of t particles. A central tool in the study of 
many-body problems is the T-matrix, describing the for- 
ward scattering between the impurity and a spin-^ atom 
at total 2D momentum q and energy e. In vacuum, the 
T-matrix is simply related to the scattering amplitude 



'7o(q, e) = /oo(e - eq/2)/m. 
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FIG. 2: The first two terms of the geometric series for the 
propagator of the impurity dressed by a particle-hole excita- 
tion. 



with the free particle dispersion Cq = /2m. The pres- 
ence of the Fermi sea couples the center of mass and 
relative harmonic oscillator modes and the T-matrix de- 
pends on harmonic oscillator quantum numbers ni, n2 
and n'^, n'2 of incoming and outgoing particles. The for- 
malism for the T-matrix in the presence of the Fermi sea 
was recently developed in Ref. [19] where it was demon- 
strated that the full T matrix, T^/",^, can be written in 
terms of a T-matrix depending only on center of mass 
quantum numbers A^, N' . For details see Ref. [19] and 
the Supplementary material. 

In the limit of weak attractive interactions, the impu- 
rity exists as a quasi-particle dressed by a cloud of major- 
ity atoms. This impurity problem has been widely stud- 
ied theoretically in both three [12, 20] and two dimen- 
sions [14-16, 21, 22] using variational and diagrammatic 
methods. In 3D, the accuracy of this analytic approach 
has been validated by both exact quantum Monte- Carlo 
calculations [23], and also by precise experiments with 
quantum simulators employing optically trapped ultra- 
cold quantum gases [24]. Here we provide an extension 
of the diagrammatic approach to the experimentally rel- 
evant quasi-2D setting. 

Consider the propagation of the impurity, initially in 
the state of harmonic oscillator quantum number n^. 
Interactions with a particle from the majority Fermi 
sea may change the state into a final state with quan- 
tum number Uf. Conservation of parity further restricts 
fii — rif to be even. In the single particle- hole approxima- 
tion, the propagator of the impurity then takes the form 
of a matrix equation 



G(p,e)= [Go-i(p,e)-S(p,e) 



(7) 



obtained by summing the geometric series illustrated in 
Fig. 2. Note that the diagrams are formally identical 
to those considered in 3D [20], however each propaga- 
tor is now assigned a harmonic oscillator quantum num- 
ber which is summed over intermediate states. The bare 
propagator is given by the diagonal matrix Gq^ (p, e) = 
^nn' / — Cpn + ^0] and the self energy is 

^n,n, (P, e) = ^ ^n^n (P + ^ + C^n) UF^iq, n), (8) 



where Cpn = + nuj^- The Fermi function np takes 
the value 1 if the state with momentum q and harmonic 
oscillator quantum number n is occupied in the t Fermi 
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FIG. 3: (color online). Energy ep / ep of an impurity interact- 
ing repulsively with a Fermi sea as a function of interaction 
strength. The energy of a quasi-2D impurity is shown for 
^f/oJz — 0.1 (solid line). Additionahy we show the 2D results 
of Ref. [16] using the often used assumption a2D = 1/y^mes 
(dotted), and also using Eq. (1) as definition of a2D (dashed). 



FIG. 4: (color online). Energy of the attractive polaron (solid 
line) and molecule (dashed) in the 2D regime and in quasi-2D 
[ep/^Jz = 1/10). The molecule energy has been shifted by 
to reflect that a particle has been removed from the Fermi 
sea. The inset shows the quasi- 2D energies of the main figure 
with the binding energy subtracted. 



sea, and zero otherwise. The energy of the polaron corre- 
sponds to a pole of ^(p, e). From Eq. (7) the energies of 
both the attractive and the repulsive polaron may then 
be obtained. 

The energy of the repulsive polaron as a function of 
interaction strength is shown in Fig. 3. As opposed 
to the attractive polaron, our calculation demonstrates 
that the energy of the repulsive polaron depends solely 
on the interaction parameter kFCi2D- This is to be ex- 
pected as \ep\ <^ oJz and thus the energy should match 
previous studies of the repulsive polaron in 2D [15, 16]. 
Furthermore, Fig. 3 demonstrates that in order to com- 
pare the predictions of these 2D theories to experiments 
with quasi-2D Fermi gases, one should not match the 
binding energy of the 2D theory to that of the quasi-2D 
system. Instead the 2D scattering length a2D defined in 
(1) and calculated by using the low energy quasi-2D scat- 
tering amplitude in Eqs. (2)- (5) should be used. Note 
that while the energy of the repulsive polaron depends 
solely on kFCi2D^ the decay rate depends sensitively on 
the energy difference between the repulsive polaron and 
the lower lying states [25] and is thus non-universal. 

As the strength of interactions is increased, the im- 
purity binds a particle from the Fermi sea to form a 
molecule. As in the case of the polaron state, the 
molecule is in turn dressed by particle-hole pairs. The 
molecule dressed by one particle-hole pair has been stud- 
ied in 3D [26-28] and in 2D [14] by variational and dia- 
grammatic methods. Here we extend the diagrammatic 
method of Combescot et al [26] to include harmonic con- 
finement. To this end, we note that the confinement 
does not change the structure of the diagrams needed 
for the molecule energy. The difference in the present 



problem is that all fermion propagators are assigned a 
harmonic oscillator quantum number and the T-matrix 
depends on these. Assuming ep < oj^ such that in the 
non-interacting Fermi gas only the ground state of the 
harmonic oscillator potential is occupied, the energy of 
the dressed molecule is obtained when the equation 

^„.„.n3 ^ ^Tn;n2(q_k, e + e^o " ek„3)[l-nFT(k, ng)] 

E [1 - --t(k', <)] - V E ^ 

+ V E E [1 - --T(k', n',)] (9) 

n^'n^ ^k q'k' ^k' J 

has a solution. The vertex H includes all diagrams oc- 
curring in atom-dimer scattering in a quasi-2D geometry 
[29]. The sum on is up to kp while the sum on is 
over all possible momenta. We use the notations = 
-e+ek+^^2 and£^JJj^/q = -e-eq+ek+e(^+ek+k'-q+^^2. 
For further details see the Supplementary material. 

In Fig. 4 we display the energy of the attractive po- 
laron and the molecule in the 2D limit and in quasi-2D. 
As may be seen, the polaron and molecule energies cross 
at a very small angle and thus the position of the polaron- 
molecule transition depends sensitively on the confine- 
ment. Computing the position of the polaron- molecule 
transition as a function of the strength of confinement, 
^f/^z^ yields the phase diagram in Fig. 1. As the con- 
finement is weakened from the 2D limit (or the density 
increased), the transition is seen to move towards smaller 
€b/^f- In particular, in the experiment of Ref. [10] with 
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FIG. 5: (color online). The energy of the attractive polaron 
in quasi-2D on the 3D Feshbach resonance. The dotted line 
is the 3D limit ep = -0.6066/it [12, 20]. 



parameters ujz = 27r x 78.5 kHz and = 27r x 10 kHz 
the transition was found at cb/^f ^ 6.8 while our theory 
predicts the transition to occur at cb/^f ^ 6.2 (see Fig. 
1). This should be contrasted with the 2D prediction of 
cb/^f ^ 10 at the polaron molecule transition [14]. 

Our formalism is also useful for studying the dimen- 
sional crossover from 2D to 3D in polarized Fermi gases. 
To illustrate this, we consider the energy of the attractive 
polaron on the 3D Feshbach resonance, where the bind- 
ing energy of the weakly bound dimer takes the universal 
value cb = 0.2446^^ [18]. Our results are shown in Fig. 
5. Spin-^ atoms in the n-th excited harmonic oscillator 
level become available whenever the chemical potential 
/i^ > nujz' The polaron energy is seen to develop a cusp 
each time a new state becomes available due to the in- 
crease in the density of t states. The polaron energy 
is seen to rapidly approach the 3D limit of an impurity 
dressed by one particle-hole pair. 

We remark on the resemblance between the above ap- 
proach and the study of the ground state of a spin po- 
larized 3D Fermi gas near a narrow Feshbach resonance 
[24, 30, 31]. The narrowness is due to a large effective 
range which enters the problem in addition to the scat- 
tering length and the Fermi wave vector. For fixed scat- 
tering length, an increase in the effective range leads to 
weaker interactions and a smaller binding energy. In the 
quasi-2D problem, the length-scale Iz plays the role of an 
effective range: Indeed, similarly to our results in Fig. 
1 for increased Z^, Refs. [24, 30, 31] found the polaron- 
molecule transition to be shifted towards smaller binding 
energies for increased effective range in 3D. 

To conclude, we have provided a theory for the highly 
imbalanced Fermi gas in quasi-2D. Experimentally, the 
phase diagram Fig. 1 may be mapped out by varying the 
strength of the harmonic confinement uj^ or the density 
in an ultra-cold quasi-2D Fermi gas. Upon entering the 
molecular phase, the quasi-particle residue of the impu- 



rity abruptly drops to zero [27], which leaves a clear sig- 
nature in the line-shape of RF-spectra [10, 13, 24]. Sim- 
ilarly, the energy of impurities in ultra-cold Fermi gases 
have been extracted in quasi-2D from ARPES spectra, 
obtained via momentum resolved RF-spectra [10]. 

Finally, we point out that while we considered the ex- 
perimentally more relevant quasi-2D problem, our results 
are also of relevance to the analogous problem of quasi- ID 
Fermi gases [32-35]. For example, similar to our results 
for the quasi-2D Fermi gas, one expects that the phase di- 
agram of the ID imbalanced gas provides an inadequate 
description at stronger coupling. 
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Supplementary material 

(Dated: March 1, 2012) 

Jesper Levinsen and Stefan K. Baur 
Cavendish Laboratory, J J Thomson Avenue, Cambridge, CBS OHE, United Kingdom 

T-matrix in the medium 

Consider the scattering of two particles at 2D momentum q and energy e. The T-matrix describing forward 
scattering in the medium in the presence of harmonic confinement was recently calculated in Ref. [19], we include 
it here for completeness: For incoming and outgoing harmonic oscillator quantum numbers ni, n2 and 77,2, the 
T-matrix describing forward scattering in the medium takes the form 

r:i:^(q'e) = V2^Z. J2 C^Z'C"N''iKMMiO)'T'N,N'{ci,e), (10) 

NN'riru'^ 

in terms of relative and center of mass quantum numbers n^, and N' . The Clebsch-Gordan coefficients for the 
change of basis are given by [36] 

r— j- min(nr.,n2) . \ / \ 

m=max(nr-— ni ,0) 

and the harmonic oscillator wavefunction takes the value 

r (n-l)!! 

at the origin. Tnn' is renormalized by the use of the two-body T-matrix as [19] 

Tnn' (q, e) = 7^"^ (e - Nuj^ - e^/2)5NN' - T>nn' (q, e) , (13) 
with the renormalized polarization operator given by 

T, (r.A- \^ r^^i^2 renins n^t (Q " k, m ) + (k, 722 ) 

t fcq— k fck ^zyf'l ~r il2) 
nin2k ^ ^ ^ 

and expansion coefficients Un^m = (— l)^^^^^/^(n— l)!!(m— l)!!/\/n!mI for n, m even non- negative integers, otherwise. 
The Fermi function ni?^(k, ni) is 1 for a particle in the t Fermi sea, otherwise. For the polarized gas considered in 
this Letter, the Fermi function ni?^(k, 712) vanishes. 

Equation for the dressed molecule 

We now construct the sum of diagrams needed to obtain the energy of the molecular state dressed by one particle- 
hole pair in the presence of harmonic confinement. Formally, the diagrams are identical to those considered in Ref. 
[26] for the 3D problem, although in the present problem all fermion lines have additional harmonic oscillator quantum 
numbers. For simplicity, we restrict ourselves to a molecule at rest. We also require ep < (j^z such that hole propagation 
always proceeds in the lowest harmonic oscillator state. Our results are easily generalizable to finite momentum and 
occupation of higher bands. 

The bound state energy corresponds to a divergence of a two-particle vertex. As argued in Ref. [26], the energy 
of the molecule dressed by one particle-hole pair necessarily lies below the energy of the bare molecule. Thus it 
suffices to consider the vertex H describing the sum of diagrams containing a maximum of two forward propagating 
t atoms as well as the impurity atom. The integral equation satisfied by H is depicted in Fig. 6: In the vertex H 
the initial interaction is between the impurity and the t atom with quantum number ni. After this interaction, this 
particle either closes its own loop (the last two diagrams on the r.h.s. of Fig. 6) or participates in another vertex 
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FIG. 6: (color online). The integral equation (9) satisfied by the vertex H. Circles indicate that the interaction inside the 
vertex is initially between the two marked particles. The loop on the last two diagrams yields a summation over q^ 



H. Taking into account only the homogenous terms (any pole in the non-homogenous term lies above the dressed 
molecule energy) finally results in the integral equation (9) in the main text. 

Additionally, we have verified that a variational approach with the trial wave functions for the polaron 



l^) = EMo„|FS,iVt) 



nn' m 
kq 



<'"clq-k„4n'kCt™q|FS,iVt) 



(15) 



and the molecule 



(16) 



kk'q 



yields the equations (7) for the polaron and (9) for the molecule quoted in the main text. Here c^j^^ denotes the 
creation operator for a particle of spin <j with in-plane momentum k in the n-th transverse harmonic oscillator mode. 
|FS,7V'^) is the non-interacting ground-state of spin-^ fermions. 



